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RESUMEN 

Se usa una estrecha corrclacion cntrc 3 propicdades de la emisi on 7 de los Es- 



tallido s de Rayos Gamma (ERGs) con corrimicnto al rojo z conocido ijFirmani et al 
l2006al) para constrenir parametros cosmologicos (PCs) en el diagrama de Hubble 
(DH) con una muestra de 19 ERGs en el amplio rango de z = 0.17 — 4.5. El 
problema de la circularidad se resuelve con un enfoque bayesiano. Encontramos 
que la cosmologia de concordancia ACDM es consistente con los datos de los ERGs 



a nivel de varias pruebas. Si suponemos el modelo A, entonces rim=0.31_QQg y 

h0.20 
-0.30 



fiA= 0.80to;3o (Ict); el caso piano csta dentro del la. Suponiendo planitud, obten- 



emos r2m=0.29^Q'Qg, y fijando Q^=0.28 obtenemos la ecuacion de estado de la 
energia oscura w ~ — 1.07^Q3g, estando cl caso ACDM {w = —1) dentro del la. 
Dado el bajo mimcro de ERGs utiles no se puede aiin constrciiir bicn la cvolucion 
de w = w{z), pcro encontramos que cl caso w{z) = —1 (ACDM) cs consistente al 
68.3% CL con los ERGs. Demostramos como un amplio rango de z's en la muestra 
usada (como cs el caso de los ERGs) mejora la determinacion de los PCs en el DH. 

ABSTRACT 

Recently, a tight correlation among three quantities that characterize the 
prompt em ission of long Gamma -Ray Bursts (GRBs) with known redshift z, was 



discovered (jFirmani et al.ll2006al ). We use this correlation to construct the Hubble 
diagram (HD) with a sample of 19 GRBs in the broad range of z = 0.17 — 4.5, 
and carry out a full statistical analysis to constrain cosmological parameters (CPs). 
To optimally solve the problem of circularity, a Bayesian approach is applied. The 
main result is that the concordance ACDM cosmology is fully consistent with the 



find f7ni=0.3llQQg and 0.801q3q (Icr); the flat-geometry case is within la. 

Assuming flatness, we find nn,=0.29tom^ and fixing f2m=0.28, we obtain a dark 



GRB data at the level of several tests. If we assume the A cosmology, then we 

hO.09 1 ri n on+0.20 

-0.08 and 17a= 0.80_o 30 

hO. 

-0.06' 

energy equation of state parameter w = — 1.07to;38, ^-^^ ACDM model {w = —1) 
is within la. Given the low number of usable GRBs we cannot yet constrain well 
the possible evolution of w = w{z). However, the case w{z) = — 1 (ACDM) is 
consistent at the 68.3% CL with GRBs. It is shown also how a broad range of z's 
in the used sample improves the determination of CPs from the HD, which is the 
case of GRBs as distance indicators. 

Key Words: COSMOLOGICAL PARAMETERS COSMOLOGY: OB- 
SERVATIONS DISTANCE SCALE GAMMA RAYS: 
BURSTS 

1. INTRODUCTION for our understanding of the universe and its fate, 

rp, . , , . , 4-- 1 1 mainly those related to the nature and physics of 
1 he impetuous advance m observational cosmoi- , , , r ii n 1 
r,, 11, the dark energy (hereafter DE) responsible for the 
ogy of the last decade has prompted new challenges ,1 . r i • 
current accelerated expansion of the universe. Stim- 
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of high energy physics, the unification of gravity and 
quantum physics, gravity beyond Einstein relativity, 
and extra dimensions. At the same time, new astro- 
nomical measurements to constrain DE parameters 
are being developed with the crucial goal of improv- 
ing quality and reducing syste matic uncertainties 
due to astrophysical effects (e.g. iLinder fc Huterer 
20051) . The main task for the new observational stud- 
ies is to tell us whether DE can be interpreted in 
terms of either a cosmological constant A (the mini- 
mal case) or something more complex and changing 
with time, such as scalar fields. In this endeavor, 
alternative and complementary methods and experi- 
ments are mandatory in order to increase the feasibil- 
ity and rigor of the results. The use of long gamma- 
ray bursts (GRBs) as cosmological distance indica- 
tors is gaining popularity as a promising method for 
constraining the cosmological parameters related to 
the dynamics of the universe. Here we present new 
advances on this method. 

As the most powerful explosions in the uni- 
verse, long GRBs are of great interest for observa- 
tional cosmology because they can be detected up to 
very high redshifts, the current record with spectro- 
scopic determinatio n being GRB 050904 at z = 6.29 



( Kawai et al. I l2006h . Ghirlanda, Ghisellini, & Laz- 



zati (2004a) have discovered a tight correlation be- 
tween the rest frame coUimation corrected energy 
and the peak energy E'pk of the i/Fi^ prompt emission 
spectrum for a sample of GRBs with known z. The 
use of this correlation has proved to be very useful 
as a method for "standardizing" the GRB energet- 
ics and its further application for constructing the 
Hubble diagram. 

The "Ghirlanda" relation has been already used 
to obtain cosmological constraints, after applying ad- 
equate approaches to overcome the problem that, 
due to the lack of a local (cosmology-independent) 
calibration, this relation actually depends on the cos- 
mological parameters that we pretend to constrain 
(Ghirlanda et al. 2004b; Firmani et al. 2005; Xu, 
Dai, & Liang 2005; Ghirlanda et al. 2006). As a 
result, the accelerated expansion of the universe at 
the present epoch was confirmed independently with 
GRBs. Interestingly enough, the marginal inconsis- 
tency of the "gold set" of Typc-Ia supcrnovae (SNIa 
hereafter) with the simple flat-geometry Friedmann- 
Lemaitre-Robertson- Walker (FLRW) cosmology in- 
cluding the cosmological constant (A-cosmology) 
(e.g. Riess et al. 2004; Alam, Sahni, & Starobin- 
sky 2004; Choudhury & Padmanabhan 2004; Jas- 
sal, Bagla, & Padmanabhan 2005; Nesseris & 



Perivolaropoul os 2005b) is eliminated when the GRB 
data are added ( Firmani et al.ll2005 : Ghirlanda et al 
20061 ). 



It is important to stress that GRBs (i) are de- 
tected from redshifts much higher than SNIa, and 
(ii) some degeneracies in determining the cosmo- 
logical parameters are reduced if the observational 
sample displays a broad r ange in redshif ts, attain- 



ing high values of z (e.g. IWeller fc Albrccht 20 ^ 



Linder fc Hutereij 20031: Nesseris fc Perivolaropoulo; 



2005aHGhisellini et al.ll2005f ). In ?4 this question will 
be amply discussed, showing by concrete examples 
why a sample broad in redshifts improves the deter- 
mination of the cosmological parameters. 

The E.y-Ep]f_ relation takes into account the 
GRB coUimation-correctcd energy, £'.^=£'iso(l- 
cos^j), where Eiso is the isotropic-cquivalcnt energy 
and 9j is the semi-aperture jet angle. The determi- 
nation of this angle is model dependent. For the 
uniform jet model in the standard fireball scenario, 
6j can be determined by the time tbrcakwhen the af- 
terglow light-curve becomes steeper. For a homo- 
geneous circumburst medium 9j oc ij^^^j^i^ (e.g. Sari, 
Piran, fc Halpern 1999), while for a wind circum- 

burst density pro file decreasing as 



^1/4 
J °^ ''break 



(|Nava et al.l l2006l ). Note that to estimate the jet 



angle from tbreak one must also assume a specific 
value of the density of the circumburst material, and 
the efficiency to convert the fireball kinetic energy 
into the radiation emitted during the prompt phase. 
Liang fc Zhang (2006; see also Firmani et al. 2006a) 
found a purely empirical multi-variable correlation 
among Eiso, i?pk and tbreak (which is then model- 
independent and assumption- free). They used this 
correlati on to constrain the cosmological parameters 
(see al solXulbOOStlGhirlanda et al.ll2006f) . 
In 



Firmani et al 



(j2006al . hereafter Paper I) 
we have searched for empirical correlations among 
7— ray prompt quantities alone. In the GRB rest 
frame the considered quantities were the bolomet- 
ric corrected Liso and i?iso, the spectral peak energy 
i?pk, and the light-c urve variability V an d duration 
To. 45 (as defined in iReichart et al. 2001 . for more 
details see § 2 below). In Paper I a ACDM cos- 
mology with r2m=0.3, riA=0.7, h=0.7 was assumed 
to calculate luminosity distances. For the sample of 
19 GRBs, for which all the above quantities can be 
defined, we have found a very tight multi-variable 
correlation among three quantities, namely Liso 
£:iif22--a49^ Within the framework of the firebaU 
scenario, the tightness of the correlation is explained 
by its scalar nature. We have also es timated the cor- 
relation among -Biso, B-pk and tbreak ( Liang fc Zhang! 
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20061) for the 15 GRBs with measured ibroak of our 
sample, and have proved that the Liso-^'pk~7o.45 re- 
lation is as tight as the i?iso-^'pk-ibreak one. 

Similarly to the "Ghirlanda" (or "Liang & 
Zhang") correlation, the iiso~^'pk~'7o.45 correlation 
can be used as a cosmic ruler for cosmographic pur- 
poses. From a practical point of view, the great ad- 
vantage of the iiso^-E'pk^To 45 correlation is that it 
involves quantities related only to the j—ray prompt 
emission. Thus, the establishment of this correlation 
avoids the need to monitor the afterglow light-curve 
in order to derive tbrcak which enters both in the 
"Ghirlanda" and in the "Liang & Zhang" correla- 
tion. 

In this paper we analyze in detail the cosmo- 
graphic application of the Liso-^'pk~2o.45 relation 
by using the current dataset. The Hubble diagram 
is constructed up to redshifts as high as z = 4.5. 
We also describe the Bayesian formalism to solve the 
'circularity problem' and compare it with other for- 
malisms. Note that this problem, at least formally, 
is also present for SNIa sa mples, as is the cas e for the 
recent SN Legacy Survey ( Astier et al. 2006h . Thus, 
the Bayesian formalism can also be used to obtain 
improved cosmological constraints from SNIa sam- 
ples. 

The GRB sample and the iiso~^'pk~'7o.45 corre- 
lation are presented in § 2. The changes of the cor- 
relation with cosmology are analyzed in § 3, where 
we test the robustness of such a correlation for cos- 
mographic purposes. In § 4 we present our approach 
to parametrize the evolution of DE, and we discuss 
the degeneracies present in the set of dynamical cos- 
mological parameters We also discuss the Bayesian 
formalism for solving the circularity problem, com- 
paring it with the conventional approach. In § 5 
we present the constraints on the parameters that 
describe the geometry and dynamical evolution of 
the Universe obtained with the sample of 19 GRBs. 
The summary and a brief discussion on the current 
shortcomings and the future of the method presented 
here are given in § 6. 

2. THE SAMPLE AND THE Liso-£^pk-To.45 
RELATION 

The sample of GRBs with known redshifts and 
with the necessary observational information avail- 
able was presented in Paper I. The rest frame 
LisQ-iJpk^To 45 correlation presented in Paper I in- 
volves: 

• the bolometric corrected isotropic energy Eiso, 
computed in the rest frame 1 — 10^ keV energy 
range; 



• the peak energy Epy^ of the i^F^ prompt emission 
time integrated spectrum; 

• the time T0.45 spanned by the brightest 45% 
of the total light curve counts above the back- 
ground and calculated in the 50-300 keV rest 
frame energy rang^. 

In addition to the spectroscopically measured 
redshift z, the observational data required to esti- 
mate Liso, fi'pk and To. 45 are the peak flux P, the 
fluence F, the spectral parameters of a given spectral 
model (in most cases the Band et al. 1993 model) 
and the light curve (to estimate To*!!)- The uncer- 
tainties in these observables are appropriately propa- 
gated to the composite quantities Liso, i'pk a-nd T0.45 
under the assumption of no correlation among the 
measured errors. Note that all the above quantities 
(except z) are obtained exclusively from the 7-ray 
prompt emission of the burst. 

In Paper I we have used a flat-geometry A cosmo- 
logical model with r2m=G.3, Ha—O-T, and h=0.7 to 
calculate the GRB luminosity distances, cJl, and to 
estimate Liso- Then, for the 19 GRBs with available 
observational data a multi-variable regression anal- 
ysis, taking into account errors in all the variables, 
provided the following best fit: 



52.11±0.03 



10 



f To.- 



pk 



1.62±0.08 



V lo"- 



102-37keV, 
-0.49±0.07 

erg s" 



(1) 



For a detailed discussion of this correlation, the er- 
ror estimates, the comparison with other correlations 
and its interpretation we refer the reader to Paper I. 

3. A COSMOLOGICAL TEST FOR THE 
-/-iso-£^pk-To.45 CORRELATION 

A preliminary cosmological test concerns the sen- 
sitivity of the iiso^-E'pk-ro.45 correlation to the dy- 
namical cosmological parameters. In Paper I we 
have assumed the currently conventional cosmolog- 
ical model. Now, we will analyze how the correla- 
tion and its scatter change from one cosmology to 
another. For each A FLRW cosmology character- 
ized by (rJiii! ^^a) we perform the multiple variable 
regression analysis on the dataset, using the same 
method described in Paper I. In this way, the (best 
fit) iiso~£'pk~7o.45 correlation, its relative scatter 
and the corresponding Xr value for each cosmology 



•^We used the recipe proposed bv lReichart et ah 1 1 I2OOII) to 
transform the observed energy range to the rest frame, and 
the time binning of HETE-II, 164-ms (see Paper I). 
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Fig. 1. Contours at 68.3%, 95.5%, and 99.7% CL's 
obtained by projecting on tlie (fim, ^^a) plane the 
Liso-i?pk^T'o.45 relation Xr derived from the fit of the 
GRB data at each value of the (fim, ^^a) pair. The star 
shows where the Xr reaches its minimum, while the cross 
indicates the concordance cosmology. This plot shows 
that the relation I/iso-^?pk-ro.45 is sensitive to cosmol- 
ogy, so that it may be used to discriminate cosmological 
parameters if an optimal method to circumvent the cir- 
cularity problem is used. The diagonal line corresponds 
to the flat geometry cosmology, the upper curve is the 
loitering limit between Big Bang and No Big Bang mod- 
els, and the lower curve indicates the division between 
accelerating and non-accelerating universes. 



are obtained and can be used to assign a prob ability 
to the (r^in, riA)-pair (jChirlanda et al.|[2004bl ). 

In Figure [1] we show the resulting contours at 
68.3%, 95.5%, and 99.7% confidence levels (here- 
afther CL's), which measure how the Xr (related to 
the scatter) of the iiso-i?pk-7o.45 relation changes 
with the cosmological parameters. Figure [1] reveals 
an important sensitivity of the scatter on cosmol- 
ogy and shows the rather surprising result that the 



Fig. 2. Contours of constant m in the (Sim, ^a) plane 
m is the power of iJpk in the Liso-i?pk-Tb.45 relation: 
Li. 

Figure 1 



oc E^Tq "g . The other curves in the plot are as in 



smallest scatter is found for (T^m, ^a) = (0.31, 
0.78), close to the concordance model (Tim, ^k) = 
(0.28, 0.72) which falls deep inside the 68.3% confi- 
dence level region. This simple and direct ('scatter- 
scanning') formalism for constraining cosmological 
parameters does not optimize the use of the avail- 
able information and is par ticularly sensitive t o 
the loitering line singularity ( Firmani et al.|[2005 ). 
However, it already shows the potentiality of the 
Liso-Ep^-ToAb relation for cosmographic purposes. 
This encourages us to use a more sophisticated for- 



malism in order to obtain more accurate cosmologi- 
cal constraints (see § 4.3). 

Figures [2] and [3] illustrate, respectively, how 
the powers m and n of the -/jiso^£'pk^'7o.45 relation 
(Liso oc i?™To7"5) change, in the (fim, ^^a) plane. 
The lines in Figures [2] and [3] arc not to be confused 
with CL contours on the cosmological parameters. 
Notice the behavior of the isocontours near the loi- 
tering curve, where the dependence of cJl on the cos- 
mological parameters becomes singular. The expo- 
nents of the i^iso~£'pk^2o.45 relation do not change 
dramatically in a wide range of (0^, Q.^) values, 
even if these changes are significantly larger than the 
small standard deviations of the exponents obtained 
in the fits (see e.g. Eq.[T]). We hope that these results 
can help for the theoretical interpretation of the ob- 
tained correlation, indicating the (rather small) range 
of the allowed m and n values. 

4. CONSTRAINING THE COSMOLOGICAL 
PARAMETERS IN THE HUBBLE DIAGRAM 

4.1. Cosmological Models with Dark Energy 

The accelerated expansion of the universe is often 
explained by the dominance in the present-day uni- 
verse of a self- repulsive medium (DE) with an equa- 
tion of state parameter w = Pde/pdec^ < —1/3. 
The simplest interpretation of DE is the homoge- 
neous and inert cosmological constant A, with w = 
— 1 and (Ode = PA =const. The combinations of 
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Fig. 3. Contours of constant n in the (fim, ^iC) plane, 
n is the power of 70.45 in the Liao-Spit-To.45 relation 
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cx E^tJTq 2^ ■ The other curves in the plot are as in 



different cosmological measurements tend to favor 
models where DE is A and (f^^, »a) «^(0.28, 0.72) 



(the s o-called concordance model, e.g. ISoergel et al, 
20031 : iTegmark et al-lbOOilSeljak et al.ll2005h . Ne^ 



ertheless, it is important to note that, due to a vari- 
ety of degeneracies in the parameter space, there are 
not yet any reliable joint constraints to the complete 
set of cosmological parameters, even after combining 
different cosmolog ical probes and data samples (e.g. 



Bridle et al.ll2003f ). Different probes can even lead 
to constraints which are not in complete agreement 
among them (when treated separately) , as is the case 
of WMAP observations of the CMB and the "gold 
set" of Type-la SNe (Jassal et al. 2005). Note that a 
more recent analysis based on the SN Legacy Survey 
has reduced this apparent discrepancy by favoring. 



the s imple flat A cosmological model (lAstier et al, 
20061) . 

Through samples of "standard candle" objects, 
such as Type-la SNe or GRBs, it is possible to 
construct the Hubble diagram and, by comparing 
the data points with the model curves (for different 
choices of Jim and I^a), to constrain these cosmolog- 
ical parameters. It is clear that, allowing w to have 
values different from —1, or even evolving, increases 
the number of free parameters to fit. Up to now, 
the existing datasets do not allow to fit together all 
these parameters. The most common approach is to 
fit only a couple of cosmological parameters, keep- 



ing all the others fixed. Such an exercise is in any 
case important since, for instance, the cosmological 
constant explanation of DE fa ces serious theoretical 
problems (s ee for reviews e.g. PadmanabhanI 2003t 
Sahnill2004[ ). Therefore, alternative scenarios, where 
w is different from —1 or even variable with z, have 
been proposed and extensively investigated. 

According to the approach mentioned above, we 
proceed here in three stages. First, we constrain the 
two parameters, (f2m, Ha), of the (minimal) A cos- 
mology {w = —1), and further check whether the 
concordance model (implying flat geometry) is sta- 
tistically consistent with the constraints. Next, we 
generalize to models with w =const (static DE), but 
assuming a flat geometry in order to have only 2 fit- 
ting parameters, flm and w =const. Finally, we gen- 
eralize to evolving (dynamical) DE models, where w 
changes with z according to a parametric form, as- 
suming a flat geometry and fim = 0.28. In the two 
last stages, with some redundancies, we again check 
whether the concordance model is statistically con- 
sistent with the constraints, i.e. whether w = —1 = 
const is within the 68.3% CL region. For the dy- 
namical DE models, we explore also how much the 
observational constraints favor the case of an evolv- 
ing or a static w. Note that any parametrization of 
w{z) is limited and arbitrary. 

To model an evolving DE we use a rather general 
parametrization for w proposed by Rapetti, Allen, 
& Weller (2005): 



w{z) = Wo + wi- 



Zt 



(2) 



where the parameter wq gives the present-day (i.e. at 
z = 0) equation of state; wi = Woo — wo gives the in- 
crement of w from the present value to z = c» and zt 
is a redshift transition scale. Note that zt should not 
be confused with the transition redshift where the 
expansion goes from decelerating to accelerating). 
The derivative of w{z) at present is Wq = wi/zf. 
The evolution of the Hubble parameter is given by 



= Hi [H^{i + zf + n^fiz) + nk{i + zf 

where = 1 — Jim — ^de, 

/(2) = (l + 2)3(i+^«-)e-3(^° 

and 



9 



I - at 
l-2at 



In 



1 - at 



a(l — 204) + at 



(3) 
(4) 

(5) 



where a is the scale factor and at ~ 1/(1 + z;) is the 
corresponding transition scale factor. 
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Fig. 4. Dependence of the equation of state parameter 
w on z as described by Eq. For the plot, wq — — ! 

is assumed. Values of wi=2, 1, and —1 (from top 
to bottom respectively) and of zt=0.5, 1.0, and 1.5 (see 
labels inside the panel) for each Wi are used. The dotted 
lines are the tangent lines of each curve at 2: = and 
represent the linear approximation. 



The simple linear approxi mation commonly 



lear approxi mation commonly 
used in previous works (e.g. iRiess et aL 2004 
iFirmani et al. 2005 ). is obtained by making zt arbi- 
trarily large and assigning a given value to w'q , which 
in this case is the slope of the w(z) function. The 
parametrization of iLinded (|2003h is recovered by set- 
ting Zt = I {at = 1/2). Figure 4 shows the family of 
parametric curves given by Eq. ([2]). Here wq = —1, 
wi = —1, 0, 1 and 2 and zt = 0.5 (short-dashed), 1.0 
(long-dashcd) and 1.5 (point-dashed). The dotted 
lines show the linear approximation at present. 

Three aspects related to the task of constraining 
w{z) are worth of mention. 

1. Methods based on the construction of the Hub- 
ble diagram with a given class of standard can- 
dles provide the primary source of information 
on the evolution of DE Q which is expected to 
become dominant only at low rcdshifts (< 1). 

2. If the redshift range of the sources is small, in 
particular limited to low z's, then the DE pa- 



rameters and evolution can be constrained only 
in a limited way (see § 4.2). 

The constraints on w{z) depend on the (arbi- 
trary) assumed parametrization for w{z). In 
fact, the space of all possible parametrizations is 
infinite-dimensional. By choosing "reasonable" 
parametrizations, both in the physical sense and 
in that of the limitations of the observational 
data, the main information we may intend to 
derive refers only to general aspects such as 
whether there is evidence or not of DE evolu- 
tion, and what is the direction of this evolu- 
tion (e.g. iLinder fc Hute rer 2005!) • Adequate 
parametrizations are those with a minimum 
number of parameters but allowing the widest 
range of variation of w over the z range in which 
w'{z) is best con strained by the given class of 
stand ard candles (jUpadhve. Ishak fc Steinhardt 
2OO5I ) . For the parametrization of Eq. ([2]) , the 



*We notice that besides of the methods based on the lu- 
minosity distance-2 diagram are also the methods based on 
the angular diameter distance— 2 diagram (e.g. the baryonic 
acoustic oscillations). 



smaller is zt, the larger the allowed change for 
w{z) at low z's, where observational data are 
available (see Figure S]). 

4.2. The Hubble Diagram for High Redshift Objects 

We now discuss some aspects related to the Hub- 
ble diagram used to constrain the cosmological pa- 
rameters. GRBs are the natural objects for extend- 
ing cosmographic studies up to very high redshifts, 
and thus for inferring the behavior of DE, in particu- 
lar, whether and how it evolves. A fundamental issue 
is, therefore, to understand all the power of the infor- 
mation which can be extracted from using the GRBs 
as standard candles extending up to very high red- 
shifts. In particular, one should be aware of the sev- 
eral degeneracies (correlations) that appear among 
the cosmological parameters at different redfshifts. 

To study such degeneracies and to understand 
the shape of the CL's in the parameter space, it is 
instructive to explore the behavior of the luminosity 
distance c?l at different redshifts z in a given cosmo- 
logical parameter space. Consider first the A cosmol- 
ogy, where the parameters are and Q\. In coun- 
terclockwise rotation, the stripes shown in Figure [5] 
represent the regions of the (f^nu ^^a) plane where 
c^L varies by ± 1% for z = 0.5, 1, 1.5, and 3, respec- 
tively, assuming that each stripe passes through the 
fiducial point (Om, f2A)=(0.33, 0.77) (sec below the 
reasons for this choice). 

The stripes in Figure [5] show that the degeneracy 
(correlation) between and ^Ia varies with z. This 
has immediate implications for cosmographic meth- 
ods based on luminosity distance measurements. 
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Fig. 5. Regions of ±1% variation around lines of con- 
stant cJl in the (fim, ^a) plane, assuming that each line 
passes through the fiducial point (f7m,SlA) — (0.33, 0.77). 
In counterclockwise rotation, the regions are at redshifts 
0.5, 1.0, 1.5, 2.0, and 3.0, respectively [in the electronic 
version of this paper z — 0.5 (magenta), 1.0 (red), 1.5 
(green), 2.0 (cyan), and 3.0 (blue)]. This plot illustrates 
the degeneracy between the parameters Slni and JIa (A 
cosmology) and how this degeneracy does change with 
z. The ellipses are contours at 68.3%, 95.5%, and 99.7% 
CL's for the fit to a A cosmology from the GRB Hubble 
diagram, using a Liso-ii'pk-To.45 relation supposed to be 
known, and therefore fixed and cosmology-independent 
(see text for more details). Note that the main orienta- 
tion of the ellipses is along the "stripes" with z ~ 1.5, 
which corresponds roughly to the typical redshifts of the 
GRB sample. The other curves in the plot are as in 
Figure 1. 



Taking into account the measurement uncertainties, 
a specific "standard candle" determines a range of 
luminosity distances c^l and consequently a stripe 
on the {flia, JIa) diagram. For a sample of stan- 
dard candles characterized by a small range in red- 
shifts, the corresponding CL's in the {flm, ^a) dia- 
gram will be very elongated (high degeneracy) and 
will have the major axis oriented in the direction 
of the stripe of the average redshift of the sample. 
Therefore, a counterclockwise rotation of the CL's is 
expected when the average redshift of the standard 
candle sample used to derive the CL's increases. This 
easily explains why the CL's derived by using SNIa 
data are elongated and oriented approximately along 
the direction of the z ^ 0.6 stripe (see Figure [5]), 
while the contours derived using our GRB sample. 



of larger average redshift (z ^ 1.5) are more "verti- 
cal". Note that, although our GRB sample contains 
a factor of 10 fewer objects than SNIa, it produces a 
comparatively narrow contour region, thanks to the 
broad distribution of redshifts of the GRBs in the 
sample. 

Figure also shows that the width of the stripes 
(i.e. the uncertainty in {Qrm ^a) associated to a 
given luminosity distance) decreases for larger z's. 
This is a consequence of the topology of the surfaces 
of constant c^l : at low redshift the surface is a gently 
tilted plane, at high redshifts the surface is more 
warped, and there appears a "mountain" with a peak 
close to ilni^ 0.0 and JIa^ 1. As a consequence, the 
stripes at high z's are curved, and at very high z's 
they surround the "mountain peak". Note that, as a 
consequence of the increasing slope of the c^l surface, 
the width of the stripes at high redshifts becomes 
narrower for large OA-values. 

From Figure [5] we conclude that in order to re- 
duce the degeneracy and improve the accuracy of 
the constraints of and Oa by using the luminos- 
ity distance method, the sample of observed sources 
should span a range of redshifts as large as possible. 

The fiducial point, («,„, Oa) = (0.33, 0.77), and 
the CL contours in Figure O have been calculated in 
the following way. We began by using arbitrary trial 
values for and Oa to define a fiducial (unique) 
Liso-£^pk~7o.45 relation. Further, we calculated the 
X^'s in the whole (Om, ^^a) plane using such fidu- 
cial Liso^Ep]^-TQ,4Q relation to assign the luminous 
distance to each GRB of known z. If the minimum 
of the x^'s was smaller than the Xr corresponding 
to the trial (Om, ^a) values, then the (Om, ^^a) 
values corresponding to the minimum Xr were used 
to define a new fiducial iiso-i?pk~2o.45 relation in a 
new iterative step. This procedure was repeated un- 
til convergence. The CL's correspond to the 68.3%, 
95.5%, and 99.7% {la, 2a and 3a) probabihties pro- 
vided by the x^ statistics. The procedure should be 
considered here as a (naive) simulation to constrain 
Oin and Qa by using a unique and well calibrated 
Liso-£'pk-2o.45 relation. Interestingly enough, the 
convergence (Om, ^^a) values in our exercise lie close 
to those of the concordance cosmology. However, we 
remark that this procedure is based on incorrect as- 
sumptions; it is introduced here only for heuristic 
reasons. 

With an analysis similar to that applied in Fig- 
ure 1, we also study the behavior of c^l in the dia- 
grams (Om, u'=const) and (wq, wi) for flat geome- 
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Fig. 6. Same as Figure [5] but in the (f7m, w) plane 
for a flat cosmology with static DE. The fiducial point 
is (fim, w) = (0.45, —2.00). Taking the vertical regions 
of the stripes, the redshifts are 0.5, 1.0, 1.5, 2.0, and 3.0 
from right to left, respectively [in the electronic version of 
this paper z = 0.5 (magenta), 1.0 (red), 1.5 (green), 2.0 
(cyan), and 3.0 (blue)]. For low values of flm, w is almost 
independent from Qrn, while the opposite happens for 
high values of Qui- The dependence of the Q^-w degen- 
eracy on z is weak. As in Figure O the bent ellipses are 
the contours of CL from the corresponding GRB Hubble 
diagram, using a Liso--Epk-T'o.45 relation supposed to be 
known, and therefore fixed and cosmology-independent. 

try cosmological models. In the latter case, we have 
further assumed that flm ~ 0.28 and used Eq. [5] 
with zt = 0.5 and 1.5 to describe the evolving DE. 
Figures [H] and [7] show the regions of (iL=const for 
z = 0.5, 1, 1.5, 2, and 3 (see details in the figure cap- 
tions) assuming in each case that the center of each 
stripe passes through a given fiducial point. 

The fiducial point in each case is [(r^m, w = 
0.45,-2.00), (wo, wi= -1.00,1.08) for zt=0.5, and 
[wo, wi= -1.01, 1.61) for zt=1.5]. The CL's in the 
(ilm, w) and (ujQ, wi) plane, represented in Fig- 
ure ini and Figure [7] respectively, were computed fol- 
lowing the same procedure described above for the 
(ilm, ^^a) plane. The la, 2<t, and 3a CL's in Fig- 
ure [S] and [7| were provided by the corresponding 
statistics. Figures [S] and [7] show the degeneracies be- 
tween w =const and flm, and between wq and wi (for 
zt=0.5 and 1.5), and how these degeneracies depend 
on z. 

To summarize: the study of the di^{z) surfaces in 
the different planes helps us to understand the ori- 




Fig. 7. Same as Figure [5] but in the (lOo, wi) plane for 
a flat geometry cosmology with dynamic DE and f2m = 
0.28. The evolving w{z) is parametrized according to 
Eq. ([2]) with 2:t=0.5 (upper panel) and zt — 1.5 (lower 
panel). The fiducial points are {wo, Wi) = (—1.00, 1.08) 
and (wo, Wi) — (—1.01, 0.61), respectively. Taking the 
lower part of the stripes, the redshifts are 0.5, 1.0. 1.5, 
2.0 and 3.0 from left to right, respectively [in the elec- 
tronic version of this paper z — 0.5 (magenta), 1.0 (red), 
1.5 (green), 2.0 (cyan), and 3.0 (blue)]. As in Figure [5] 
the ellipses are the contours of CL from the correspond- 
ing GRB Hubble diag ram, using a LiBo~Ep]^—To_4^ re- 
lation supposed to be known, and therefore fixed and 
cosmology-independent . 



entations of the CL regions for different samples of 
cosmological probes, characterized by different aver- 
age redshifts. This study makes intuitively clear the 
need to have probes distributed in a large range of 
redshifts. This in turns implies that SNIa and GRBs 
complement each other in a natural way. 
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4.3. The Bayesian Formalism 

Now wc will explore how the correlation 
iiso-£'pk-To.45 can be used to constrain cosmological 
parameters through the Hubble diagram. In the pre- 
vious section we introduced the concept of a unique 
well calibrated Liso--E'pk-T'o.45 relation; however, this 
is not the present case. In fact the iiso-£'pk-?o.45 
depends on the assumed cosmology. Therefore the 
crucial issue in this undertaking is what has been 
called the "circularity problem" : we attempt to con- 
strain the cosmological parameters using a correla- 
tion which is cosmology-dependent. This problem 
arises because, due to the lack of detected low— z 
GBRs, the £iso-£'pk-'7o.45 correlation can not be 
calibrated at low redshifts, where the flux is not 
affected by a specific cosmology. Another way to 
calibrate this kind of correlation is with a sample 
of high-red shift GRBs in a considerably small red- 
shift bin. Ghirlanda et al have calculated 
that ~ 12 GRBs with z € (0.9, 1.1) can be used to 
calibrate the "Ghirlanda" relation with a precision 
higher than 1%. This number might be reached in a 
few years of observations mainly due to the fact that 
the jet break time measurement (which enters in the 
"Ghirlanda" correlation) requires a time-consuming 
follow up campaign of the GRB optical/NIR after- 
glow. We estimate that a similar number of GRBs 
might also be used to calibrate the Liso-Ep]^-ToA5 
correlation. Fortunately enough, as the latter cor- 
relation only relies on prompt emission information, 
we should expect to collect few tens of GRBs with 
a low redshift dispersion in a few months, provided 
that an adequate 7-ray instrument acquires the rele- 
vant prompt emission information, namely the light 
curve and a broad band spectrum. 

While waiting for a sample of calibrators, ade- 
quate statistical approaches should be used in or- 
der to optimally recover cosmographic information 
from the cosmology-dependent points in the Hub- 
ble diagram. The B ayesian formalism presented in 
Firmani et al. ( 20051 ) is currently the most suitable 
method for this purpose and we apply it here for 
constraining cosmological parameters by using the 
Liso-i?pk~To.45 relation. The basic idea of such for- 
malism is to find the best-fitted correlation on each 
point fl of the explored cosmological parameter space 
[for instance fl = and to estimate, us- 

ing such a correlation, the scatter x^i^^^) on the 
Hubble diagram for any given cosmology fl. The 
conditional probability P{^l\Cl), inferred from the 
X^(r2, Q) statistics, provides the probability for each 
U given a possible Jl-defined correlation. By defining 
P'(f2) as an arbitrary probability for each f2-defined 



correlation, the total probability of each fl, using the 
Baycs formalism, is given by 



p{n) = / p{n\n)p'{n)dn, 



(6) 



where the integral is extended over the available Cl 
space. Note that from the observations one obtains 
a correlation for each cosmology. Therefore, P'{i^) 
is actually the probability of the given cosmology. 
Consequently such probability is obtained by putting 
P'{n) = P{n) and solving the integral Eq. ([6]). It 
should be noted that in the conventional use of the 
Bayes approach, P' is handled as a given prior prob- 
ability. Here, instead, P' and P are just the same 
probability which is solution of Eq. ([6]). 

An elegant Monte Carlo approach allows us to 
solve Eq. ([6]), i.e. to find the probability P{fl) from 
this integral equation. Wc start by determining 
the empirical correlation for an arbitrary cosmology 
Hq. The r^o-defined correlation is used to calculate 
on the Hubble diagram the probability distribution 
Po{^\flo) exp(— x^(f2, fio)/2)- From this proba- 
bility, we randomly draw the cosmology f2i, which 
is used to determine again the empirical correlation. 
Then, with the f2i-defined correlation, we calculate a 
new probability distribution pi{^) that is averaged 
with po(Cl). The result is the probability distribu- 
tion Pl{^l). From this probability, a cosmology 1^2 is 
randomly drawn and is used to determine again the 
empirical correlation. Applying this correlation in 
the Hubble diagram gives a new probability distri- 
bution p2 (f^) • The new probability p2(f^) is averaged 
with the previous ones, po{i^) and pi(ri), giving the 
probability distribution P2{^)- The cycle is repeated 
until convergence, i.e. until Pi(fi) did not change 
significantly with respect to Pi_i(fi). The conver- 
gence should be fast if the empirical correlation is 
not too sensitive to cosmology. On the contrary, 
if the correlation is strongly dependent on cosmol- 
ogy, then convergence cannot be attained with this 
method. Numerical experiments show that for the 
Liso~i^pk"To.45 relation, convergence is attained af- 
ter a few thousands of cycles with a result that is 
independent of the choice of the initial cosmology 
fig. By introducing some numerical techniques, the 
convergence can be attained after hundreds of cycles. 

The described formalism is very different from 
assuming that the correlation is known (unique and 
well calibrated) as done in the previous section. It 
is also different from scanning directly the 
rameter for all points in the plane by minimiz- 
ing the scatter of the data points around a corre- 
lation that is found in the very same Q point, and 
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therefore changes from point to point, as we did in 
§ 3. Besides, the Bayesian formalism is less affected 
than the direct method by possible discontinuities, 
like the loitering line in the (r2m,r2A) plane, and it 
avoids spurious divergences. With the Bayesian for- 
malism, the observational information, provided by 
correlations like the iiso~^pk^7o.45 one, is optimally 
extracted for cosmographic purposes as long as this 
kind of correlations remain uncalibrated at low z. 

It should be remarked that there is no formal and 
rigorous mathematical method for solving the "cir- 
cularity problem" . However, a comparison of Fig- 
ures 1 and 8 clearly shows how the constraints im- 
prove when one or another formalism is used. Xu et 
al. (2005) have also shown the much better perfor- 
mance of the Bayesian formalism as compared to the 
other methods. 

Finally, we mention th at the Bayesian appro ach 
has been also used by us (|Firmani et al. I l2006b[ ) to 
constrain cosmological par ameters from the S uper- 
nova Legacy Survey (SNLS. lAstier et al-lliooi ). The 
SNLS data are given in such a way that they also 
require a cosmology-dependent calibration. We have 
found constraints and CL contours in the (f2m,f^A) 
plane similar to those of Astier et al. (2006), who 
used the direct minimization method; if anything, 
our CL contours were slightly narrower. This shows 
that our method is working for what it has been 
designed, namely to optimize the data from quasi- 
standard candles in the Hubble diagram. 



a 1 




Fig. 8. Constraints on the {Qm, Qa) plane for a A cosmol- 
ogy from the GRB Hubble diagram using our Bayesian 
method to circumvent the circularity problem (thick-line 
ellipses) and from the "gold set" SNIa Hubble diagram 
(thin-line ellipses). The ellipses are contours at 68.3%, 
95.5%, and 99.7% CL's. The other curves in the plot are 
as in Figure 1. 



distinct trends that are not shared by the "gold" 
set (see also Nesseris fc Perivolaropouloi 2005lj : 
Firmani et al.ll2006bl ). 



5. RESULTS 

In this section we present the results on the cos- 
mological constraints obtained with the Bayesian 
formalism (§ 4.3) applied to the tight Liso~i?pk"T'o.45 
correlation defined with the 19 long GRBs dis- 
tributed in redshift to up to 4.5. Following § 4.1, we 
proceed to constrain only 2 parameters each time. 
In all the models we fix /i = 0.71. It should be em- 
phasized that our results represent a first attempt, 
still using a dataset with few numbers and with a 
quality of the observational information not yet at 
the level of SNIa. However, these results allow us to 
quantify the potentiality of the GRB iiso~^'pk~2o.45 
correlation as an independent cosmological tool. 

For comparison purposes, we will also show 
the cosmolo gical constraints p rovided by the SNIa 
"gold set" (|Riess et al.l l2004l z < 1.67). The 
latter were derived by using the standard (di- 
rect) x^"fitting procedure. It is worth to men- 
tion that the results on cosmological constraints 
recently presented by the "Supernova Legac y Sur- 
vey" group (z < 1.01; lAstier et all |2006[) show 



5.1. A Cosmology 

In Figure [8] we show the Icr, 2a, and 3(7 CL's 
(thick lines) for the f2m and Ha parameters. No- 
tice how the CL's improve with respect to those 
obtained with the simplest direct minimization 
method used in § 3 (Figure [T|). The best-fit cosmol- 
ogy (see the star symbol in Figure [5]) corresponds to 



17^=0.311^:^^, f7A=0.80Io;3^ (Ict uncertainty). This 
result is very close to the flat geometry case. The 
concordance model is well within the la CL. If the 
flat geometry case is assumed (i.e. r2tot= 1), our 
statistical analysis constrains flm = 0.29to'o6- 

The constraints on the A cosmology parame- 
ters that we obtain with GRBs alone are consistent 
with those obta ined through several other cosmologi- 
cal p robes (e.g. I Hawkins et al JEooil ISchuecker et al 
20031 ). In turn, this result gives us confldence that 
GRBs can be used as cosmological probes. 

In Figure [5] are also shown the best-fit values 
(star symbol) and CL regions (t hin lines) that w e 
obtain with the SNIa "gold set" (|Riess et al.ll2004l) . 
As these and other authors (e.g. Choudhury & 
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Fig. 9. Constraints on the (f2m, w) plane for a flat cos- 
mology with static DE from the GRB Hubble diagram 
using the Bayesian formalism to solve optimally the cir- 
cularity problem (thick-line ellipses) and from the "gold 
set" SNIa Hubble diagram (thin- line ellipses). The bent 
ellipses are contours at 68.3%, 95.5%, and 99.7% CL's. 

Padmanabhan 2004; Jassal et al. 2005; Nesseris & 
Perivolaropoulos 2005b) have shown, the "gold set" 
provides constraints on and ^Ia that are only 
marginally consistent with the concordance model 
or the WMAP GBR constraints. 

5.2. Flat Cosmology with Static (w = const) DE 

Figure [H] shows the la, 2cr, and 3cr GL regions 
on the (rim, w) plane for models with static DE and 
flat geometry, using the GRB sample (thick lines) 
and the SNIa "gold set" (thin lines). 

The degeneracy here is relevant and higher than 
the corresponding degeneracy seen in Figure [S] This 
feature is consistent with the discussion of Figures [5] 
and [6] of the previous section and has to do with the 
small rotation of constant di^ lines for different red- 
shifts. The reduction of such a degeneracy will be 
possible by reducing the GRB observational uncer- 
tainty a s well as by increasing the number of objects 
(see e.g. iGhirlanda et alll2006l . for a simulation). 

The A case {w = —1) is consistent at the 68.3% 
CL with the GRB constraints, for values of ili„ = 
0.29^0 07. The concordance model is well inside 
the 68.3% CL. For a prior r2m=0.28, we obtain 
w = — l.OT^Q'gg. Note that the A case is not con- 
sistent with the SNIa "gold set" at the 68.3% CL. 




Fig. 10. Constraints on the {woyWi) plane for a flat cos- 
mology with dynamic DE and nin=0.28 from the GRB 
Hubble diagram using the Bayesian formalism to solve 
optimally the circularity problem (thick-line ellipses) and 
from the "gold set" SNIa Hubble diagram (thin-line el- 
lipses). Upper and lower panels are for zt=Q.^ and 1.5, 
respectively. The ellipses are contours at 68.3%, 95.5%, 
and 99.7% CL's. The diagonal dot-dashed line is the 
upper limit in the (iuo,wi)-plane allowed by CMB con- 
straints. 



5.3. Flat Cosmology with Dynamical DE 

Formal constraints on (u>o, wi) (see Eq. [5]), as- 
suming flat-geometry cosmologies (and = 0-28) 
with dynamical DE, are presented in Figure [TOl Up- 
per and lower panels refer to zt = 0.5 and Zt = 1.5, 
respectively. The thick and the thin line ellipses are 
the 1(7, 2a and 3a GL regions for the GRB and SNIa 
"gold" samples, respectively. The A case {wo = — l 
and wi=0, which reduces to the concordance model 
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because of the assumption that rini=0.28), is within 
the la CL. 

The typical redshift of our GRB sample is z«1.5. 
The best information provided by GRBs on the value 
of w{z) is expected at the same redshift. Our anal- 
ysis for zt = 0.5 gives w(1.5) = — O.SI^'q, and for 
zt = 1.5 it gives u;(1.5) = -0.51?;^^. The (still) large 
uncertainties in the data and the small number of 
objects do not allow to constrain zt as a third pa- 
rameter. 

Again, note that the constraints provided by the 
SNIa "gold set" are not consistent with the concor- 
dance model values of wq and wi at the 68.3% CL. 

In general, the "gold" SNIa constraints tend to 
favor low values of wq and large values of wi, im- 
plying (i) a strong evolution of w{z) in the range 
-Z ^ 1, and (ii) a significant probability for 'w{z) 
to cross the w = —1 line (phantom divide line, see 
also e.g. Riess et al. 2004; Alam et al. 2004; Nesseris 
& Perivolaropoulos 2005b). As reviewed by the last 
authors, if observations show a significant probabil- 
ity for w{z) to cross the phantom divide line, then all 
minimally coupled single scalar field models would 
be ruled out as DE candidates, leaving only models 
based on extended gravity theories and combinations 
of multiple fields. It is therefore a key observational 
task to determine whether w{z) crosses the w = —1 
line or not. The SNIa "gold set" rejects models that 
avoid the phantom dividing line at the la CL, while 
our results with GRBs allow these models (includ- 
ing the concordance one) at the la CL, though the 
uncertainties for the latter are still much larger than 
for the former. 

Finally, we should emphasize that Eq. ^ is just 
a mathematical paramctrization for the evolution of 
w, but not a physical model of DE. Although Eq. ([2]) 
describes the evolution of w up to any arbitrary large 
z once its parameters arc determined, the changes in 
w with z suggested by the observational constraints 
are formally valid only within the redshift range of 
the observational data. For example, the constraints 
shown in Figure [10] cannot be used to extrapolate the 
behavior of w{z) as given by Eq. ^ to z's higher 
than ~ 3, and ^ 1 for the GRB and SNIa data, 
respectively. 

In fact, at high redshifts there are several ob- 
servational limits to the values of the parameters 
of Eq. ([2|). The most important is related to 
the CMB anisotropics. The CMB data require 



DE 



< 



0.1 at the redshift of recombination. 



1100 (jCaldweU fc Doranll2004[ ). For the "Rapetti" 
paramctrization that we are using (Eq. [2]) and as- 
suming flat geometry, this condition implies that 



Wo + Q.SGwi < —0.095, which is close to the general 
upper limit of Woo = wi + wq < fou nd in the anal- 
ysis of WMAP and other data sets by lRapetti et al. 
(|2005f ). The dashed line in Figure fTOl corresponds 



to this limit. Interestingly enough, the best-fitting 
point from the GRB sample in the {wq, wi) plane 
obeys the CMB constraint for zt = 0.5, being slightly 
out of this constraint for zt = 1.5. Instead, for the 
SNIa "gold set" the best-fitting points in both cases 
are far away from the CMB constraint. 

6. SUMMARY AND DISCUSSION 

Firmani et al. (Paper I) found a very tight cor- 
relation among three GRB quantities in their rest 
frame, Liso, E^y^ and T0.45. These quantities were 
calculated from the 7— ray prompt emission spectra 
and light curve, without the addition of any quantity 
derived from the afterglow, apart from the redshift. 

Here we have used this tight correlation to "stan- 
dardize" the energetics of the currently available 
sample of 19 GRBs, and to construct an observa- 
tional Hubble diagram up to the record redshift of 
z = 4.5 and independent from SNIa. Based on the 
behavior of the luminosity distance as a function of 
different cosmological parameters (§ 4.2), we have 
pointed out that samples of standard candles dis- 
tributed over a wide redshift range arc strongly de- 
sired for breaking the degeneracy of the cosmolog- 
ical parameters. To overcome the circularity prob- 
lem that arises due to the lack of a local cosmology- 
independent calibration of the Liso^Ep]i-ToA5 rela- 
tion, w e have applied a Bay esian formalism devel- 

M) and further discussed 



(200£ 



oped in iFirmani et al.l 
here. The main results on the cosmological con- 
straints are: 

• The iiso-^^pk-7o.45 correlation is sensitive to 
the cosmological parameters of the A cosmol- 
ogy (§ 3), having a minimum Xr (^^m, ^^a) 
= (0.31, 0.78), very close to the concordance 
model (Figure 1). 

• For the A cosmology, using the Bayesian for- 
malism, the best-fitting values for fim and Q\ 



are 0.31Io;o8 and 0.67Io;3o {la uncertainty), re- 
spectively. This result is very close to the flat 
geometry (Figure [8]). The ACDM concordance 
model (f2ni=0.28 and ^Ia= 0.72) is well within 
the 68.3% CL. If one assumes flat geometry, 



then we find fin 



0.29 



+0.08 
^0.06- 



For constant w models (static DE) with flat ge- 
ometry, the A case {w = —1) is consistent at 
the 68.3% CL for values of = 0.29^^;^?. 
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The ACDM concordance model is still within 
the 68.3% CL. 

• For models with dynamical DE, we have 
parametrized w{z) according to Eq. ([2]) and 
used zt = 0.5 and zt = 1.5. Assuming a flat ge- 
ometry and ri,n=0.28, the A case {wo=—l and 
wi=0, which also in this case corresponds to the 
concordance model) is again within the 68.3% 
CL. Interestingly enough, the constraint that 
the CMB data (z = 1100) provide on w{z) as 
given by Eq. ^ {wq+OMwi < -0.095), is con- 
sistent with the constraints found with GRBs. 

We conclude that the different constraints pro- 
vided by the GRB sample are consistent at the 68.3% 
CL with the ACDM concordance model. This is not 
the case of the SNIa "gold set" . Also, the GRB con- 
straints for flat-geometry models, with DE equation 
of state parameter either constant or varying with 
z, are consistent with the constant w{z) case at the 
68.3% CL, while the "gold set" SNIe are not. These 
results show that the GRB method presented here 
offers already a competitive and reliable way to dis- 
criminate cosmological parameters. 

The use of the correlation Liso-^'pk-7b.45 among 
prompt 7-ray quantities has proved to be a promis- 
ing, model-independent and assumption-free method 
for constructing the observational Hubble diagram 
up to high redshifts. The accuracy that this corre- 
lation provides in constraining cosmological parame- 
ters with the current available set of useful GRBs is 
better than that found with other correlations (either 
the "Liang & Zhang" correlation or the "Ghirlanda" 
correlation) . Most importantly, the advantage of the 
Liso~Ep\c~ToA5 correlation is that it does not involve 
any quantity related to the afterglow. 

Compared to SNIa, the GRB cosmological con- 
straints arc less accurate. This is due to the still low 
number of GRBs having the required data as well as 
to the relatively large uncertainties associated with 
these data. However, GRBs provide valuable com- 
plementary cosmographic information, in particular 
due to the fact that GRBs span a much wider red- 
shift distribution than SNIa. As discussed in § 4.2, 
some degeneracies appear when constraining the cos- 
mological parameters with samples of "standard can- 
dles" limited only to low redshifts. The results pre- 
sented in this paper are a clear proof of the poten- 
tiality of using the GRB iiso-^^pk-'7o.45 relation for 
cosmographic purposes. After the completion of this 
paper, a paper by Schaefer (2007), where a com- 
bination of several (noisy) empirical correlations of 
GRBs was used to construct the Hubble diagram up 



to z = 6.4, appeared posted in the arXiv preprint 
database service. The constraints on the cosmologi- 
cal parameters obtained in that work are similar to 
ours, though the data and methodology are very dif- 
ferent from the ones presented here. 

It is worth to mention that as more data of higher 
quality appear, some assumptions made concerning 
the cosmographic use of the iiso~£'pk-To.45 relations 
will either be accepted or refused. For example, in 
order that this relation be useful for cosmography, 
it should not evolve, or the way in which it changes 
with z should be known. It is also important to im- 
prove the quality of the data in order to reduce the 
scatter, as well as to increase the number of usable 
GRBs. So far, the best physical justification of the 
iiso-£'pk~'7o.45 relation derives from its scalar na- 
ture, which explains its reduced scatter because it is 
independent of the relativistic factor F (Paper I). A 
full physical interpretation of this relation is highly 
desirable, in particular to avoid any uncertainty con- 
cerning observational selection effect^. However, as 
an empirical relation used like a distance indicator 
tool, the main concern is related to reducing the ob- 
servational scatter. Just recall the case of the fa- 
mous Tully-Fisher relation for disk galaxies, which 
has been used as distance indicator for more than 
twenty years, though until recently its physical foun- 
dation was not clear. 

Another potential problem for high-redshift 
GRBs cosmological tool is gravitational lens- 

ing which systematically brightens distant ob- 
jects through the magnification bias and increases 
the dispersion of dist ance measurements (e.p 
IPorciani fc Madaull200ll lOguri fc Takahashil [2006 



Recently, Qguri fc Takahashil ( 2006 ) simulated the 
gravitational lensing effects on the Hubble diagram 
constructed with Swift-li ke GRBs following a reason- 
ablc luminosity function (jFirmani et al lliool . They 
showed that lensing bias is not drastic enough to 
change constraints on dark energy and its evolution. 
However, they emphasized that the amount of the 
bias is quite sensitive to the shape of the GRB lumi- 
nosity function. Thus, an accurate measurement of 
the luminosity function is important in order to re- 
move the effect of gravitational lensing and to obtain 
unbiased Hubble diagram. 

We finish by emphasizing that the ideal strat- 
egy to follow in the future is to combine the SNIa 
and GRB data sets, and to adopt the same meth- 



^In a recent paper, appeared after the first submission of 
the present one, Thompson, Rees & Meszaros (2006) have 
suggested some interesting hints to understand the origin of 
the iiso^-Bpk— 70.45 relation. 
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ods of handling these data sets of "standard can- 
dles" in order to construct their joint Hubble dia- 
gram, and thus constrain the cosmological parame- 
ters. Of course, the dominant information will be 
that of SNIa (they outnumber GRBs and the un- 
certainties on their luminosity are smaller than for 
GRBs), but GRBs provide valuable information at 
high redshifts which helps to partially overcome pa- 
rameter degeneracies and biases. This pro gram is 



carried out elsewhere (jFirmani et al.ll2006b[ ) 
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